Bistability of optically-induced nuclear spin orientation in quantum dots 
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We demonstrate that bistability of the nuclear spin polarization in optically pumped semicon- 
ductor quantum dots is a general phenomenon possible in dots with a wide range of parameters. 
In experiment, this bistability manifests itself via the hysteresis behavior of the electron Zeeman 
splitting as a function of either pump power or external magnetic field. In addition, our theory 
predicts that the nuclear polarization can strongly influence the charge dynamics in the dot leading 
to bistability in the average dot charge. 

PACS numbers: 73.21. La, 72.25Fe 
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The hyperfine interaction in solids between the elec- 
tron and nuclear spins [3| leads to the Overhauser energy 
shift S of the electron spin states, produced by the mag- 
netic dipole moments of orientated nuclear spins, often 
described in terms of the effective nuclear magnetic field 
Bn = 2S/geiJ-B- The hyperfine interaction is also re- 
sponsible for the transfer of spin from electrons to nuclei 
and may lead to a significant nuclear spin polarization S, 
if the system is pumped with highly polarized electrons 
[IjBBI • Recently, nuclear spin effects have been observed 
in optically pumped semiconductor quantum dots (QDs) 
@, M, 0, B H [lol- 1 111- IT^ . In these experiments circularly 
polarized light generated electron-hole pairs which, then, 
relax into the ground state of the dot, with electrons ex- 
hibiting a longer spin memory than the holes (which lose 
their spin polarization due to a stronger spin-orbit cou- 
pling Mm)- 

Recently, the nuclear spin orientation in optically 
pumped dots has been found to display a pronounced 
bistability in externally applied magnetic fields [13, [HI, 
[T^ . This appears as a threshold-like switching of the 
nuclear magnetic field 2 — 3T and a characteristic 

hysteresis behavior observed in the dependence of the nu- 
clear polarization on either the intensity of the polarized 
light IHI or external magnetic field [H, [l^ ■ 

In this Letter we propose a theory of the nuclear polar- 
ization bistability in optically pumped QDs. We study 
the dynamics of nuclear spins in a dot populated by elec- 
trons (el) and holes (h) which arrive into its ground state 
with the independent rates w (el) and li (h) and polar- 
ization degrees a (el) and a (h) (see Fig. HUa)). It has 
been recently noticed that nuclear polarization bistabil- 
ity may occur in the regime when light generates 100% 
spin-polarized excitons on the dot. Here, we demonstrate 
that bistability is a general phenomenon possible in a 
wide range of experimental conditions, including the non- 
resonant excitation conditions and in the regime when a 
dot often appears in a positively charged (trion) state. 
We also predict a new phenomenon caused by the bistable 
behavior of the nuclear spin orientation: the bistability 
of the dot average charge. 

In optically pumped dots, nuclear spins become orien- 
tated due to the spin flip-flops in which a single electron 
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FIG. 1: (a)Diagram of optical excitation of the dot. The case 
of strongly polarized electrons and weakly polarized holes is 
shown, (b) Spin-fiip-assisted el-h recombination accompanied 
by fiipping one of the nuclear spins. 



and one of the nuclei exchange spins via the hyperflne 
interaction. The process leading to spin transfer consists 
of an el-nucleus spin flip-flop accompanying the recombi- 
nation of a polarized electron with a spin ±| heavy hole 
in a QD carrying a neutral or positively charged exciton 
(see Fig. [TJb)). In such a process, the electron occupies 
the intermediate inverted-spin state on the dot virtually, 
since in a magnetic fleld a real single-electron spin-flip is 
prohibited by energy conservation and the Zeeman split- 
ting. The rate of the electron spin-flip recombination 
involving a single nucleus [l3| is. 



l«|V/(eL + h')- 



(1) 



Here u is a typical energy of the hyperfine interaction 
with a single nucleus, 7 is broadening of the electron 
energy level, and Wr is the rate at which the bright 
exciton recombines on the dot. The electron Zeeman 
splitting, modified by the Overhauser field Bjsi, is Cez = 
gelJ.B{B — Bn). The form of the Eq. ^ implies a feed- 
back due to the dependence of Ws on Bn, which is key 
to the nuclear spin bistability. 

The kinetic model describing the carrier population in 
the ground state of the dot is formulated in terms of the 
probabilities of its 16 allowed configurations based upon 
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the two electron and two hole spin states correspond- 
ing to the lowest el/h orbitals in the QD. We solve the 
rate equations for the populations of these states and for 
the nuclear orientation, and then find the steady-state 
magnitude of the nuclear spin polarization S. Here, we 
denote the probability that the dot is empty by n, and 
use (nf^) for the probabilities of the dot occupation 
by a single electron (hole), with the index n = +/— rep- 
resenting the spin state of the particle. We refer to these 
states as D,D^ and D^, respectively. The probabilities 
for the dot to be occupied with two electrons or two holes 
(states and ) are n_| and . The probabil- 
ity to find the dot in a dark exciton state is n^, and 
in a bright exciton state X'^^^ is n^^ . The probability to 
find the dot in a negative (positive) trion state labelled 
as X'^_ {X'^~) is (n+~) and, finally, n\_Z represents 
the dot in the biexciton state, X^Z- 

Below we list the balance equations for the dot popu- 
lation. The first two equations describe the probability 
of the dot occupation by a single carrier. 



i(l + fia)'wn + Wr'n'^ 
^(1 -f ^d')wn + Wriij^ 
— \_w + ^(1 — /iCT)ii;] n 



\w + i(l — ^a)w\ n 
i(l + tiS)Nw,n^_ 



(2) 



Both include "gains" due to the arrivals of an elec- 
tron/hole into the empty dot (Fig.[lja)) and the recom- 
bination of a charged bright exciton, and " losses" due to 
the arrival of an electron or a hole. The second equation 
also has a gain due to a possible spin-flip-assisted recom- 
bination from a positive trion X^ — > D~'^ in which the 
spin is transferred to a nucleus Fig. [ijb). This pro- 
cess is impossible for a negative trion since in the lowest 
orbital state the flip-flop is blocked by the presence of the 
second electron [16|. The probability for an el-h pair to 
recombine via spin-flip depends on the number of nuclei 
available, which leads to the term (1 -I- iiS)Nwsn~^~ in 
Eq. (21) , where S is the degree of nuclear polarization and 
TV is the total number of nuclei covered by the electron 
wave function {N ~ 10^ ^ 10^ in a typical InGaAs/GaAs 
dot). 

Equations describing the QD states ZJ-i and Z?^ are: 

= i ^^(1 ^ /^o')wv — wn^ ; 

= i^(l - AfCT)w)n'' - wn+-. (3) 

Kinetics of the the neutral bright and dark excitons 
Xj^ and X^^^ are described by 

n'^= ^{1 + ^ia)wni_, + ^{1 + fia)wn'^ 

-i [(1 - ^^S)Nws + (1 - ^i(T)w + {1- fid-)w] ni^; 
ri^^ = i(l — iia)wn^ + ^(1 + iiG)wn^^ + 7«,.n^Z 

-[wr + \{l~ ^i(T)w + \{l + ^ia)w]n-^'. (4) 



Both neutral bright and dark exciton populations de- 
crease when more carriers arrive onto the dot. The neu- 
tral bright exciton can also be created and removed due 
to the el-h pair recombination in the processes X^Z — > 
X~^ and X'^'-'- — *■ D, respectively. The dark exciton 
can decay due to the spin-flip-assisted recombination 
(X^ D) leading to spin transfer to nuclei [l6l|. 

Kinetics of the trions X'^_ and X'^~ are described by 

n'^_ = ^(1 4- fid-)'wn^ + i ^ (1 — i'(j)wn'^ 

-\wr^'\ij^~ pio)w\n\_\ (5) 
i(l + [ia)wn^ + i ^ (1 — h'a)wn'^ 



-[wr + \{l- nS)Nws + 5(1 - ^ia)w] n+-. 
Both trion populations change due to the recombinations 
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and arrival of a single additional 
charge (the ground states of the dot permit maximum 
four carriers). A positive trion can also recombine in the 
spin-flip-assisted process X'^Z^ X~'^, forbidden for the 
negative trions [16| . 

Finally, the biexciton state X^Z cannot contribute to 
the nuclear spin pumping as it decays without the spin- 

flip, x+z^xz,xl, 



\ [(1 — ii(j)'wn'^_ + (1 — fia)wn'l^ 



(6) 



The probabilities for the dot with a given nuclear po- 
larization to be in each of the 16 configurations are found 
using the normalization condition I — n + -|- n'^ -|- 



n;, and the 



steady state condition for Eqs. ([^HH). We formally write 
these equations in the form Mn = (1, 0, 0)"^, where 
the components of ft are the occupation numbers and M 
is a 16 X 16 matrix with elements determined by the co- 
efficients in Eqs. (I^Ml) and the normalization condition. 
The formal solutions for components of n are given by 
Ci^i/detM where Ci^i is the relevant cofactor of ill. 

A steady-state value for the nuclear polarization S (de- 
fined as S* = — /ji) can be obtained by substituting 
formal steady-state solutions of Eqs. ([2][6]) for a given S 
into the balance equation for the occupation numbers of 
spin up (/^) and down (/^) nuclei p^ . 



5 = /EE^^(l-;.5) 



2Swd. (7) 



It summarizes the processes leading to the nuclear spin 
pumping: S is increased as a result of the spin-flip- 
assisted recombination of and X'^~ and reduced 
due to a similar recombination process involving XZ and 
X^^ . Thus the balance between the populations of X'^ 
and Xt~ on one hand and XZ and X't~ on the other 




FIG. 2: The function I{S) for the situation where = w, cr = 
0.9, a — —0.2, X = 0.6 and z — 8 for three different powers: 
P = 0.0001, 0.0003, 0.0005. Stable roots correspond to the 
solutions of 1(5) = where H < 0. 

will eventually define the sign of the net nuclear polariza- 
tion pj^. However, an additional important contribution 
to the depolarization of the nuclei has to be taken into 
account. It arises from their mutual dipole-dipole in- 
teraction effectively leading to the nuclear spin diffusion 
from the dot into the bulk semiconductor [19| , described 
in our model by the rate Wd [l^l . 

To present the analysis of the above equations, we em- 
ploy the following parameters: 



B 



TDmax ' 



Wr 



P 



— . (8) 

ZWr 



Here B^"^ is defined through the Overhauser field Bn as 



Bn = B^"^S and, after introducing a = "f/gefJ^sB 
Eq. Il]) can be represented in the form [2ll |: 



max 
N ' 



a mm sum 
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FIG. 3: Evolution of nuclear polarization (5) and the average 
charging state of the dot (Q) for w = w,cj = 0.9 and d — —0.2: 
(a) as a function of power for x = 0.6 and various values of 
z, with the arrows indicating a forwards or backwards sweep. 
Although not shown in the figure, at high powers P « 1, both 
Q and S start to decrease due to the dot being dominated by 
the biexciton (for which the spin-flip process is blocked); (b) 
as a function of magnetic field for z = 8 and various power 
values. 



{x - S) + 



4" 



(9) 



The steady-state values of S determined by the feedback 
built into Eqs. ([T][2]) are given by the solutions of the 
equation / (5*) = 0, satisfying the condition ^ < (so- 
lutions with ^ > are unstable). Figure[2]demonstrates 
that for a fixed external magnetic field the number of sta- 
ble solutions for the nuclear spin polarization varies: it 
can be one or two depending on the incident power and 
other experimental parameters such as Wd,cr,a and the 
ratio w/il). At small powers only a single low value of 
S is possible. At high powers when two stable solutions 
appear, including one with a large S > x, the dot enters 
the regime of the nuclear spin bistability. This result 
strongly depends on the depolarization parameter z, de- 
fined in Eq. ([5]), so that in the following discussion we 
specify the range of z where a bistability occurs. 

The bottom parts of Figs, ^a) and (b) show the cal- 
culated evolution of the nuclear polarization in a dot 
for realistic magnitudes of the depolarization parame- 
ter z in the regime where electrons have a high degree 



of spin memory and arrive with the same rate as the 
depolarized holes. Fig. [3^a) contains a large hysteresis 
loop in the power dependence of S for a fixed magnetic 
field (here, x = 0.6), similar to those observed in Refs. 
10,11. The bistable behavior occurs for a wide range of 
the depolarization parameter z: 5 < z < 14. Experimen- 
tally, the evolution of S can be detected in polarization- 
resolved PL experiments on individual self-assembled In- 
GaAS/GaAs quantum dots, by deducing it from the mea- 
sured exciton Zeeman splitting. 

We also find that the bistability in S leads to a novel 
phenomenon: a hysteresis in the average dot charge, Q 
(see top parts of Fig.[3Ia,b)). This occurs when the elec- 
trons arriving to the dot have a high degree of spin po- 
larization, permitting their recombination with only one 
spin orientation of holes. Thus, an extra hole with the 
opposite spin is likely to remain on the dot, leading to, 
on average, a positive dot charge. The enhancement of 
the spin-flip-assisted recombination for large 5, removing 
such holes, will result in reduction of the charge. There- 
fore the hysteresis in S will be reflected as a hysteresis 
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FIG. 4: Evolution of S with P in the regime where w — 
O.lw for z = 0.16,0.4 and different polarizations of arriving 
electrons/holes, (a) a — 0.9, cr = —0.2. The inset shows the 
evolution of the charging state of the dot for z = 0.4 (with a 
very small hysteresis loop), (b) Same for a = 0.45, a = —0.1. 



in the average dot charge. A similar bistable behaviour 
in both S and Q can also be found if the external mag- 
netic field is varied [l^] at a fixed optical pump power, 
as shown in Fig. [Sljb). 

Figure 2] illustrates that the range of parameters for 
which the bistability can occur strongly depends on the 
ratio between the arrival rates of electrons and holes, w 
and w, as well as on their polarizations, a and a. In ex- 
periment the ratio w/w can be varied by applying an elec- 
tric field in a diode containing QDs in the intrinsic region 
d, H, [nl, [13] : because of a hght effective mass, electrons 



can tunnel out before relaxing to the dot ground state, 
which in effect reduces their arrival rate as compared to 
that of the holes. Fig. ID^a) shows the evolution of S{P) 
for w = O.lw. The dot is mainly in the state D'^ so 
that its average charge is Q ~ -1-1.8 and exhibits a weak 
power-dependence with a negligible hysteresis loop (see 
inset), despite a pronounced hysteresis loop in the nu- 
clear polarization. As seen from the figure, for such low 
values of w/w higher powers are required to pump a sig- 
nificant nuclear polarization, and the bistability in 5' is 
moved towards smaller values of the depolarization pa- 
rameter (0.2 ^ z < 0.5). Figure (Hb) illustrates that 
when the polarizations of both electrons and holes is re- 
duced by 50% the bistability can still be observed, but 
only for 0.1 < z < 0.2. 

To summarize, we have shown that for a wide range 
of dot parameters (including the number of nuclei, el-h 
radiative recombination time and nuclear spin diffusion 
rate) the polarization of nuclei in a non-resonantly op- 
tically pumped semiconductor quantum dot can exhibit 
a bistable behavior. Thus, we conclude that the nuclear 
spin bistability is a general phenomenon for dots pumped 
with circularly polarized light. In addition, we find that 
the nuclear spin polarization can also strongly influence 
the charge dynamics in the dot leading to the bistability 
of the average dot charge. 
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